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This study presents two methods of analysis for the 
supersonic oscillating cascade with subsonic leading edge. 

A relatively straightforward solution is developed for the 
slowly oscillating finite and infinite flat-plate cascades 
which provides simple analytical expressions for the unsteady 
pressure distributions. Comparison with other solutions 
is generally excellent. Some additional linear topics 
including resonance and a unique inflow condition are also 
treated. In addition a nonlinear method of characteristics 
solution for finite cascades is described which permits 
analysis of blade thickness effects on flutter. At this 
time, only the inlet and passage flow computations have 
been completed and are compared with available experimental 
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1. INTRODUCTION 



1.1 General. Unsteady flow in turbomachines can be broadly 
divided into two categories ; circumferential asymmetry in 
the flow through the turbomachine and flow perturbation due 
to the self-excited oscillatory motion (flutter) of the 
internal blading [1]. As shown in figure 1.1.1, taken from 
[2] , at least four types of flutter may occur in a compressor 
rotor. The two stalled-f lutter phenomena are encountered 
when compressor or fan blading is subjected to incidence 
angles near or exceeding the stall angle, and so are off- 
design phenomena. Subsonic stalled flutter has long plegued 
the engine designer [3] . Choke flutter and supersonic unstalled 
flutter are low-incidence instabilities occurring when the 
flow over the blade is, respectively, transonic along most 
of the chord, or supersonic over the outer span portions [3] . 
Supersonic unstalled flutter (hereafter referred to simply 
as supersonic flutter) , which has appeared only recently as 
a significant problem in turbomachinery, is the subject of 
this investigation. Primarily a concern in compressors and 
fans [1] , it is an instability born of the recent technologi- 
cal advances aimed at increasing the thrust-to-weight ratios 
of turbomachines: increased tip speeds, composite blades, 

lighter discs and blades, and the elimination of part-span 
shrouds [2,3]. Some insight into the nature of supersonic 
flutter is gained by a consideration of its characteristics: 




Figure 1.1.1 Compressor map showing boundaries for four 

types of flutter. 



1. Its flutter boundary imposes a high-speed limit on 
compressor operation [2] . 

2. It is a low-incidence, attached-flow instability that 
can occur at the design condition (figure 1.1.1). 

3. Flutter amplitude decreases with increasing back 
pressure on the rotor [2] . (For this reason the present 
analysis considers the most critical case of zero pressure 
rise across the rotor.) 

4. At the flutter condition all the blades are observed 
to flutter at the same frequency with a constant interblade 
phase angle [2]. 



5. Both single-degree-of-freedom torsional flutter (pitch) 
and two-degree-of-freedom torsional-flexural flutter (pitch 
and plunge) can occur [3,4]^ 

6. The flutter boundaries are significantly affected by 
airfoil shaps [4,5]* 

7. The primary flutter mechanism is an inviscid phenomenon. 
If the unsteady aerodynamic force lags the blade motion 

by a phase angle which is sufficient to ensure positive 
work over the entire cycle, then flutter can occur [3,4]. 

It is clear that reliable techniques for predicting the 
unsteady pressure distribution on a fluttering blade are 
required if one hopes to predict and avoid the onset of 
flutter. With this end in mind the cascade has been exten- 
sively employed as a theoretical abstraction of a compressor ^ 
rotor. The unsteady flow is considered in an annulus of 
differential radial height within which equispaced airfoils 
are radially disposed. Unwrapping such an annulus gives 
rise to the two-dimensional approximation of plane flow 
through a cascade of airfoils [1]. The cascade is termed 
linear if the airfoils are flat plates of zero thickness. 

AS seen in figure 1.1.2, a rotor does not require a super- 
sonic axial flow to have a supersonic velocity component 
relative to, and parallel to, the blade axis. Indeed this 



A Arroldi et al. in "Supersonic 

Flutter in Cascades Pratt & Whitney * gQ^^^^cur 

of 31 Hay 1975 show that chordwise bending can also occur, 




Figure 1.1.2 Supersonic cascade classification. 



distinction provides a convenient method of classifying 
supersonic cascades. If the axial velocity is subsonic, 
the cascade is said to have a subsonic leading edge; in the 
other case it has a supersonic leading edge [6] . The physi- 
cal significance of this distinction is clear from figure 
1.1.2. In the supersonic-leading-edge cascade each blade 
influences only its immediate neighbors. When the cascade 
has a subsonic leading edge, however, disturbances generated 
by the upstream portions of each blade affect the flow around 
each subsequent blade in the cascade; furthermore the flow- 
field at the underside of each blade is influenced by the 
wake disturbances of all preceding blades. Since the supersonic- 
leading-edge cascade has been thoroughly investigated [7,8,9, 

10], the more' difficult subsonic-leading-edge case, which is 
also the more interesting in practice, is considered in the 
present study. Another cascade classification is frequently 
employed: if each blade is assumed to be preceded and followed 

by infinitely many others, the cascade is said to be infinite; 
if there is an identifiable first blade — a necessary pre- 
requisite for many numerical solutions — the cascade is 
termed finite. 

A considerable amount of effort has been expended in 
attempting to predict — numerically and analytically — 
the unsteady pressure distribution over a blade in an oscilla- 
ting supersonic cascade with subsonic leading edge. A 
brief discussion of some of the more noteworthy advances 



follows . 



1.2 The Linear Subsonic-Leadinq-Edge Cascade . Verdon 



and McCune [11] formulated a linear boundary value problem 
for the subsonic-leading-edge configuration, solving it 
for velocity potential. This permitted calculation of the 
pressure along the entire upper surface of an arbitrary blade 
and along the lower surface from the leading edge to the 
point of impingement of the trailing-edge Mach wave from 
the preceding blade. In a more recent paper Verdon [12] 
formulated a second boundary value problem, in which pressure 
is the basic dependent variable, which permitted determination 
of the pressure along the entire lower surface as well. 
Nagashima and Whitehead [13] and Goldstein [14] obtained 
solutions by replacing the cascade blades with doublet dis- 
tributions. Kurosaka [6] introduced what he terms the 
passage approach in applying Laplace transform techniques 
to a slowly-oscillating cascade, obtaining closed-form 
analytical expressions for the pressure distributions of 
certain configurations. He recently [15] extended his solu- 
tion to include arbitrary frequencies of oscillation. In 
assuming that a space periodicity of the flow exists from 
one blade passage to the next, Kurosaka claims that he 
avoids complications arising from the breakdown of linear 
acoustic theory in the far field [4] . Platzer and Brix 
[16] and Platzer and Bell [17] employed the method of char- 
acteristics to solve directly for the perturbation veloci- 
ties and pressures everywhere in the flowfield. Their 
computer program has been thoroughly checked against existing 




solutions [11,12,13] and is extensively employed in the 
present study as a method of comparison for both the analy- 
tical and numerical solutions presented. Chadwick [5] used 
Sauer's solution [18] of the unsteady potential equation, 
accurate to a first order in frequency of oscillation, to 
determine the pressure distribution on a given blade from 
the leading edge to the point of the first reflection (point 
A in figure 1.2.1) of disturbances from the preceding blade. 
He obtained exact analytical agreement with Kurosaka [6] over 
that portion of the blade for the single-degree-of-f reedom 
pitching cascade. Chadwick's technique has the decided 
advantage of being much simpler than that of Kurosaka, easily 



yielding closed-form expressions for any cascade, and 
raising the question as to whether the method could be 
extended to the entire blade for the cases of both pitch 
and plunge. Finally, Adamczyk [19], in a yet-unpublished 
work, employs the Wiener-Hopf method to obtain the pressure 
distribution for the linear cascade with subsonic leading 
edge. 

1.3 The Nonlinear Subsonic-Leading-Edge Cascade. Attempts 
to deal analytically or numerically with an oscillating 
cascade of blades having thickness are not nearly so numerous 
as their linear counterparts. Carrier [20] provided an exact 
theory for a single oscillating wedge. VanDyke [21,22]. 
also .treated the isolated wedge as well as an oscillating 
airfoil of arbitrary profile. Teipel [23,24] employed an 
ingenious method-of-characteristics approach to numerically 
obtain the unsteady pressures on a single oscillating profile 
of arbitrary shape, obtaining good agreement with the single- 
blade solutions of Carrier and VanDyke. Teipel ' s method 
assumes potential flow with negligible entropy changes 
across leading-edge and trailing-edge shocks. Kurosaka and 
Edelfelt [4] have obtained numerical results for an isolated, 
oscillating, parabolic-arc airfoil and are presently working 
to extend their solution to a cascade of such airfoils. 
Chadwick [5] employed Teipel ' s method to obtain the unsteady 
pressure distribution over the preinterference zone of an 
arbitrary blade in a cascade of wedges. His results indicate 



that nonlinear thickness effects significantly alter the 
unsteady pressure profile of an arbitrary blade. 

1-4 The Present Investigation . The present work has two 
major parts: The first is analytical and treats the linear 

cascade while the second is numerical, dealing with a non- 
linear configuration. 

In the linear effort a closed-form analytical solution 
for the pressure-difference distribution along a blade in 
a slowly-oscillating cascade is obtained. The solution 
technique is that developed by Chadwick [5] in the pre- 
interference zone (figure 1.2.1) of a pitching cascade, in 
which Sauer's solution of the unsteady linearized potential 
equation is employ-e’d. In the present work Chadv/ick's solu- 
tion is extended through all reflection zones and into the 
blade's wake. In addition, the solution is obtained for 
two degrees of freedom (pitch and plunge) . Two standard 
configurations — Verdon's cascades A and B [12] — are con- 
sidered, one from a finite-cascade approach and the other 
as an infinite cascade in which Kurosaka's passage approach 
is employed. Extensive comparisons are made with existing 
solutions [6,12,19,17]; and several special topics, including 
cascade resonance and a unique inflow condition, are briefly 
treated. 

The nonlinear effort is an evaluation of the feasibility 
of extending Teipel's unsteady method of characteristics 
to a cascade of airfoils with arbitrary profiles. Teipel 



[23,24] represents the nonlinear flowfield by a harmonic, 
unsteady, velocity potential superimposed on a steady one, 
and discontinuities in the flowfield variables across the 
shockwaves emanating from the airfoil are determined from 
the isentropic Rankine-Hugoniot relations. These concepts, 
coupled with the characteristic grid and compatibility relations 
derived by Teipel, permit determination of the steady and 
unsteady preinterference pressures acting on the second 
blade of a cascade of parabolic-arc airfoils. The numerical 
results are compared with experimental data from tests con- 
ducted with a pitching cascade by Fleeter and Riffel [25j 
and with linear theory. An improved method for determining 
the jump conditions across the bow shock is presented and 
the computational procedure is carefully outlined to facili- 
tate extension of the technique to the entire flowfield of 
a completely general cascade. 

1.5 Some Comments on Notation . In the linear portion of 
this investigation it is the unsteady, dimensionless, 
pressure-difference distribution along a particular cascade 
blade that is sought. This pressure difference is a complex 
quantity given the symbol £^C and represents an amplitude 
function of distance along the blade. The harmonic oscilla- 
tion factor is always implicitly understood. In 

addition it is desirable to obtain an expression for 

AC at the instant that the blade in question achieves 
c> 



maximum up displacement, referred to as the initial position 

This expression is obtained by multiplying through by the 

factor e and so it is indicated by the sumbol AC e-^/" 

P 

In the nonlinear portion of this study only the local 

pressures — C upper and lower — are sought. They are 
P 

also complex quantities and it is convenient to express 
them as a magnitude and phase angle relative to the motion 



of the blade under study; the symbols 
respectively, are employed. 



and 
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2 . 



THE LINEAR LOW-FREQUENCY CASCADE 



2.1 General . A cascade of flat plates with chord c 
oscillating harmonically with frequency OJ in a steady, 
uniform flow U is said to oscillate at a dimensionless 
reduced frequency of oscillation k defined by 




( 2 . 1 . 1 ) 



If k 



2 



is small when compared with k , then the cascade is 



termed a low-frequency or slowly-oscillating cascade. Although 
low-frequency flutter is not common, it is worthy of investi- 
gation for several reasons, the major one being that it is 
the only case for which one can hope to attain a relatively 
simple, closed-form expression for the unsteady pressures 
on a given blade. Such expressions provide insight into the 
flutter phenomenon, as evidenced by the identification of 
a unique inflow condition in a subsequent chapter of this 
study. Also, the low-frequency cascade provides an important 
limiting case for any general-frequency solution, numerical 
or analytical. 

2 . 2 Problem Formulation . 

2.2.1 The Linearized Potential Equation. The flow 



variables in the cascade flowfield are assumed to be 
composed of a freestream component and a much smaller 



perturbation component 






( 2 . 2 . 1 ) 



The cascade is assumed immersed in a perfect each 

particle of which maintains constant entropy as it moves 
through the flowfield 



0^ 

D-t 



-O 



( 2 . 2 . 2 ) 



In addition the flow is assumed inviscid and irrotational , 
guaranteeing the existence of a velocity potential 

V = (2.2.3) 



The continuity and momentijm equations for the flow are, 
respective ly 



|£--vy> v:v=o 
-22. + ZP. =5 

D-t ^ 



(2.2.4) 

(2.2.5) 



Introducing the velocity potential and inteR;rating (2.2.5) 
produces, after linearization 



at a-x. A. 

Using the equation for the freestream 



O 

speed 



( 2 . 2 . 6 ) 

of sound 






(2.2.7) 

the continuity equation can be linearized to give 



^ ^ \ -V- T-J -O (2.2.8) 
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From (2.2.6) one obtains 



and 



St 



sx 



fsV 1 



(2.2.9) 



r 

Uxs-t 



•xj 



sHzi 1 

J vJ 



( 2 . 2 . 10 ) 



Substitution of (2.2.9) and (2.2.10) into (2.2.8) then 
produces the linearized unsteady potential equation in 
tv7o dimensions 



3t» 



2.±L£!^ _ H 

ixst ~ ^ ^ sx*- ^ 






(2.2.11) 



x^7here M is the freestream Mach number 



M = 



XJ 



( 2 . 2 . 12 ) 



2.2.2 The Linearized Pressure Coefficient . A generalized 
form of Bernoulli’s equation follows directly from the 
momentum equation (2.2.5) 
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iL - 



= CONST, 



(2.2.13) 



Introducing the perturbation quantities in the form of 
(2.2.1) and ignoring the products of perturbation 
quantities one obtains 



where the prime notation has again been dropped. Thus, the 
dimensionless pressure coefficient is given by 




2.2.3 The Linearized Boundary Conditions . The surface 
of a body placed in the flow can be described as a function 
of space and time 






9-L ^ 7 . 




(2.2.14) 



rcw.-tw o 



(2.2.16) 



where r is the position vector of points on the surface. 

The condition that the flow at the body surface be everywhere 
tangent to the surface can be expressed by 




Figure 2.2.1 Two-degree-of-freedom oscilla-tion. 



at the surface. If the body is a flat plate oscillating 
harmonically in plunge and/or pitch, as shown in figure 2.2.1, 
then the describing functions are respectively 

^ ^ ^ (2.2.18) 

Ki (2.2.19) 

where h and are the maximum displacements and b 

o o 

is the pitch axis location on the x axis . Applying the 
generalized flow tangency condition (2.2.17) to (2.2.18) 
and (2.2.19) yields the boundary conditions for plunge and 



pitch respecitvely 



= -ho <x-t ^ = o (2.2.20) 

-|^ = i.oocx-b^e"'^*' o.-t^rO (2.2.21) 

2.2.4 Nondimensionalization. It is convenient to intro- 
duce nondimens ional variables into the governing equations 
to simplify and generalize them. The distance variables 
(x,y,b) are normalized by chordlength c , velocity potential 
by the product Uc and time by the quantity c/U . It is 
also convenient to set 

C. 1. RA.DIKN (2. 2. 2-2) 

and to separate the velocity potential into a complex 
amplitude and phase angle 

(2.2.23) 

Thus the potential equation (2.2.11), pressure coefficient 
(2.2.15) and boundary conditions (2.2.20,21) become 

= ^ (2.2.24) 

r ;»x*- ^ 






(2.2.25) 
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(2.2.26) 



~ a-t (2.2.27) 



where , 



yg^=r M^-1 



(2.2.28) 



and C and ^ are complex amplitude functions of x and 

sr 

ik t 

y only with the factor e omitted. 



2.2.5 Cascade Geometry. The physical parameters chosen 
to describe the cascade are the vertical distance between 
the blades d and the horizontal distance between the 
leading-edge Mach waves m . Both parameters are normalized 
by chordlength. The cascade can also be described by 
solidity ^ and stagger angle 0 , which are shown in 
figure 2.2.2. These parameters are related as follows 



in this study are Verdon's cascades A and B [12] , which 
have become standards in the field and are depicted in 







(2.2.29) 






(2.2.30) 



d 



where yS is given by (2.2.28). The linear cascades employed 



figure 2.2.2. 
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2.3 The Isolated Airfoil. Consider a single flat plate 
oscillating with two degrees of freedom as shown in figure 
2.2.1, and define 

(2.3.1) 

Then s is constant along any left-running (upper surface) 
Mach wave, while s is constant along right-running (lower 
surface) waves. Sauer [18] showed that the potential 
equation (2.2.24) has the following left-running and 
right-running solutions, respectively, for a slowly- 
oscillating airfoil 

(2.3.2) 

(2.3.3) 

where the g and h functions are determined from the 
boundary conditions (2.2.26,27). Requiring that the 
velocity potential be continuous across the leading-edge 
Mach waves, one has 

^ ^ O (2.3.4) 

Enforcing flow tangency at the upper surface of the airfoil 
yields, for plunge and pitch respectively 




Figure 2.2.2 Cascade geometry. 
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huCs>= — S, 



(2.3.5) 



V.*,ts'»=s-^^( ^-.-b^*s) (2.3.6) 



And so on the upper surface 



PL; 



PI; 



= ^ s. 



(2.3.7) 



^[t (2.3.8; 



Thus the upper-surface pressure coefficient follows from 
(2.2.25) 



PL; 



PI; 



C“ = -Sife. 



(2.3.9) 



,= -^ (2.3.10) 



where the subscript indicates the single airfoil. Similarly 
for the lower surface 



Cp,= -C“ (2.3.11) 

2 . 4 Continuity of the Velocity Potential Across a Mach Wave 
Consider two reflection zones separated by the Mach wave 
s = s^ , as shown in figure 2.4.1, and let the variables z 
and z respectively measure distance tangent and normal to 
the Mach wave. The velocity potential in zone II is simply 






y 




igure 2.4.1 Variation of the velocity potential along 
the boundaries of reflection zones. 

at in zone I plus an additional disturbance 

^ (2.4.1) 

nee the tangential velocity component remains unchanged 
ross the wave, (2.4.1) implies that 9 ^^ is a function of 
alone, along the Mach wave 

4-OT So (2.4.2) 

it since z is itself constant along the Mach wave, then 
, must have the form 



CONST. +or s = 



(2.4.3) 



That the constant is zero follows from the physical con- 



straint that the velocities remain finite in the wave 
itself, and so 






(2.4.4) 



which is in agreement with the conclusion reached by 
Landahl [26] . 

2 . 5 The Finite Two-Degree-of-Freedom Cascade. As an 
elementary means of introducing the method of analysis, 
the pressure-difference distribution is found for the second 
blade of a finite version of cascade A. This analysis also 
affords the opportunity to compare the finite-cascade .results 
obtained with numerical finite-cascade data [17] and analy- 
tical results [6] of existing solutions. It is emphasized 
that the pitch results through equation (2.5.27) were first 
obtained by Chadwick [5] . The reflection zones referred to 
are pictured in figure 2.5.1, and the pressure results are 
given in terms of the isolated-blade pressure coefficients 
(2.3.9,10). Before proceeding with the analysis, some 
additional varieibles are introduced for blade two 



which are related to those of blade one by the translation 





(2.5.1) 



relations 







(2.5.2) 



s=:cr-vm ssscr-t-vw-^ 



2ySd 



(2.5.3) 



2.5.1 Zones I and lA. In zones I and lA Sauer's single- 
blade solution applies 

(2.5.4) 



T* T 

(2.5.5 



and in zone lA the boundary conditions take the form 



PL: 



CL-t ^=0 



(2.5.6) 
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a 02 
a ‘i 



= -[^1-*. LV?(or-b)] <tt. »t=iO (2.5.7) 



where yuc is the interblade phase angle by which the second 
blade leads the first. Applying the boundary condition at 
the upper surface of each blade and using that 



= Vi.co') = O (2.5.8) 



gives , 


for plunge and 


pitch 
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^,(S> = o 




(2.5.9) 


PI: 




= ^ -»r B 


(2.5.10) 


PL: 


h,l« =-V S 


h^^(a-^=: 4-C 

r 


(2.5.11) 



PI: ^ (2.5.12) 

h,^ (O') = - 4-^ Ci M Vj 

where A, B, C and D are constants of integration. By 
enforcing continuity of the velocity potential across the 
leading-edge Mach wave of blade two, the constants are 
evaluated and the velocity potentials deterinined in terms 
of the isolated blade potential 9^“(2.3.7,8) 



PL: 



i>? - i>'t 






(2.5.13) 



= 9^r 

PI; (2.5.14) 

- 1^1 ^ ^ - ^,(v«MV+y5v>,N^*^ + ji) 

where 0 is a constant that need not be evaluated because 

2 

it becomes the coefficient of a term in k when the 
pressure coefficients are formed.^ Thus, on the upper 
surface of the second blade, in terms of the isolated-blade 
pressure coefficient Cp^ (2.3.9,10) 

PL; = Cp, (2.5.15) 

PI; (2.5.16) 

r 

2.5.2 Zones IIB and VIB . In zone IIB the velocity 
potential is that of the isolated blade (2.3.3) 

(2.5.17) 

while in VIB the potential acquires an additional 



^The symbol 0 is used throughout section 2 to indicate 
a constant which need not be determined. 



TQ 



right-running disturbance 






(2.5.18) 



2.5.3 Zone II. In zone II a right-running reflection 
occurs and, using the translation relations (2. 5. 2, 3) to 
shift to the blade-two coordinate system, one has 



PL : = 9S^ V>z ■*' 

PI: 

+^,(5*-) -t k[h2(o=) ^ \ 



(2.5.19) 



(2.5.20) 



where 



^nCo'> sVijCo") =o 



(2.5.21) 



Applying the flow tangency condition at the lower surface 
of blade two and using the Sauer functions to construct the 
velocity potentials yields 



PL: 





(2.5.22) 



5=[-^»m +^elM»+ (2.5.23) 

-VO" cvn-v ^4 -V igm 



And the pressure coefficients follow from the potentials 
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(2.5.24) 



PI: 



(2.5.25: 



Thus the pressure difference can be computed through zone 
II, and by multiplying by the factor e the pressure- 
difference amplitude is obtained at the moment the second 
blade reaches maximum displacement 

PL: r^-2Cp,Cl-e”’^> (2.5.26) 
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(2.5.27) 
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2.5.4 Zone III . In zone III a left-running wake distur- 
bance is experienced, and after applying the translation 
relations 



PL: 



95“ = 5i/ + (2.5.28) 
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- S (W1 - e^/*)] (2.5.29) 
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where 



= VijCl) =. o (2.5.30) 

Since zones III and VIB are separated by the wake of the 
blade rather than the blade itself, the wake boundary con- 
ditions are used between the two zones. The first of these 
demands pressure equilibrium across the wake 

Cp^ = a.-t (2.5.31) 

which leads to the equations 

= S>T (Cp,) (2.5.32) 

(2.5.33) 

where the right-hand sides of (2.5.32,33) refer to the 
steady and frequency terms respectively of (2.3.9,10). 

The second wake condition requires continuity of the normal 
velocity component through the wake 

a-t w=ro (2.5.34) 

which leads to the equations 






(2.5.35) 



V»1 (S> -+• =0 (2.5.36) 

Equations (2.5.32,33,35,36) are solved for the zone-III 
disturbances 

PL: ^j(s')=0 Kj(0) = --^ (2.5.37) 

PI: ^jCs)=-^CS-0 (2.5.38) 

2.5.5 Zone V. In zone V an unknown disturbance is added 

-4- ^sC^> v?^h3((r) + VigC^) 

(2.5.39) 

which is determined by enforcing flow tangency at the lower 
surface of blade two in zone V 



PL: 



^gCcr) ~ O hs 



hsCo-) = - 



(2.5.40) 



PI: 



(o'— I -tVh^ 

hs(O^) = ^ 



(2.5.41) 



which leads directly to the pressure difference for the 
last portion of the blade 



PL: 






^ i 



(2.5.42) 
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(2.5.43) 



<,^ t. \ 



2.5.6 Results and Conclusions. The finite results for 
the second blade are compared with four existing solutions 
in figures 2.5.2 and 2.5.3 for three interblade phase angles: 
Kurosaka's solution [6] is a closed-form, analytical, 
infinite-cascade solution; Platzer's [16,17] is a numerical, 
method-of-characteristics solution for the finite cascade; 
and Adamczyk's [19] and Verdon's [12] solutions are numeri- 
cal, infinite-cascade results. It is interesting to note 
that the finite results agree exactly with Kurosaka's 
infinite solution over the entire blade for the plunge case 
and through zone II for the pitch case. Over the zone-V 
portion of the pitching blade — the portion exposed to the 
wake disturbances from blade one — the two solutions agree 
numerically, as seen in figure 2.5.3, but are not analytically 
the same. In fact, Kurosaka's solution is singular for 
0 (in-phase oscillation of the blades) over the last 
portion of the blade, while equation (2.5.43) behaves well 
for that case. All five solutions were found to be in 
excellent agreement for interblade phase angles larger than 
15°, but for small angles considerable disagreement is 
evident, as seen in figure 2.5.2. 
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Figure 2.5.2 Comparison of the finite-cascade 
solution with existing solutions. 
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Figure 2.5.3 Comparison of the finite-cascade 
solutions with existing solutions. 



Chadwick [5] compared his finite and infinite results 
for the pitch case over the in-flow portion of the blade and 
noted that although the finite and infinite pressure- 
<iiffsrence (ACp) distributions are the same, the local 
pressure distributions (Cp, C^) may in fact differ. The 
author found this phenomenon occuring over the remaining 
portion of the pitching blade also. Figure 2.5.4 shows that 
the local pressure distributions for the present finite 
solution and Verdon's infinite solution do not agree as 
well as the pressure-difference distributions. Some con- 
clusions can be drawn from the analysis at this point: 

1. For low frequencies of oscillation the finite cascade 
rapidly approaches the infinite-cascade solution, achieving 
good agreement with the pressure-difference distribution 

of the latter at the second blade. 

2. For small interblade phase angles, the present solution 
agrees with Kurosaka's Laplace transform solution for the 
infinite cascade but considerable differences with the other 
published solutions are noted. 

3. Although the present finite-cascade solution and 
Verdon's infinite-cascade solution for the pressure-difference 
distributions are in agreement for large interblade phase 
angles, the local pressure distributions resulting from the 
two solutions are seen to disagree. 

2 . 6 The Infinite Two-Degree-of-Freedom Cascade. In order 
to demonstrate the applicability of the above analysis to 




Figure 2.5.4- Comparison of the finite and infinite 

local pressure distributions. 



the infinite cascade, it is applied to an infinite version 
of Cascade B. Since the advantage of an identifiable first 
blade is lost, a condition specifying the space periodicity 
of the flowfields between the blades is derived. The pitch 
results through equation (2.6.17) were first obtained by 
Chadwick [5] , and the reflection zones referred to in the 
solution are depicted in figure 2.6.1. 
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2.6.1 The Flow Periodicity Condition. The periodicity 
(f the flow from blade passage to blade passage was first 
.dentified by Kurosaka [6] and is the heart of his passage 
ipproach to the infinite cascade. An alternate derivation 
ollows. 

Consider the corresponding zones I and lA, and assume 
•hat the potentials are the same there except for a complex 
actor K 



^ ( 2 . 6 . 1 ) 



lince each potential must satisfy the boundary condition of 
•ts corresponding blade 
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( 2 . 6 . 2 ) 
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it is clear that 



< = e'v^ 



(2.6.4) 



and so 



95? (>,?■) = vS?" (J,*!) e-> 



(2.6.5) 



The same result follows for any pair of corresponding zones 
on adjacent blades — for .pitch and plunge — and is termed 
the flow periodicity condition. 

2.6.2. Zones I and lA. The analysis in zones I and lA 
proceeds much as in the finite case (section 2.5.1). The 
velocity potentials have the form 




( 2 . 6 . 6 ) 




(2.6.7) 



and the boundary conditions yield 
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Vines') = ^ s ~ ^ A s -H e: 

A /S 

PL: (2.6.10) 

h,A(cr> = 

h,(S) (^ -4-y^s-v^M‘Cs.')-^ G* 

PI: (2.6.11) 

h,^tO')= - ^ W^^cr) & '-/< H-ySM^3>o^ -► H 

The constants A through D are evaluated using, respectively, 
continuity of the velocity potential across the leading-edge 
Mach wave of the (n+l)st blade and the flow periodicity 
condition 






(2.6.12) 
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The pressure coefficients, in terms of those of the isolated 
blade (2.3.9,10), follow 
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(2.6.15) 



2.6.3 Zones II through IV. Finding the additional 



reflections through zone IV is accomplished by routine use 
of the flow tangency condition (2.2.26,27), and so only the 
pressure results for the (n+l)st blade are given below 

PL; (2.6.16) 

PI; ‘ ^ (2.6.17) 

O ^ ^ <2.^J 

PL: 2^<l<^<l-vn (2.6.18) 

PI: ^ . (2.6.19) 

2.6.4 Zones V, VI, and VII. In the first three zones 
above the wake of the nth blade 



S 6^=1 ( 2 . 6 . 20 ) 

0^1 = 9^^ + ^ ^ t (2.6.21) 

«s-” = 

where the g^ and h^ Sauer functions are known 
analysis in zone IV and 



from the 



= hgCO - o 



(2.6.23) 



^7C\-vv>):s V>,0-Vn) = 0 (2.6.24) 

2.6.5 Zones VA and VIII. In the zones straddling the 
wake of the (n+l)st blade, using flow periodicity and 
results from zones IV and V 

*l (2.6.25) 



where 
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n h^(0 = O 
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The wake boundary conditions — pressure equilibrium and 
normal velocity continuity - are enforced between zones 



VA and VIII 



«.-t ^ = o 
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giving rise to two equations in three unknown Sauer 
functions 



,?5 » - -<?7 -c?8 = S'T [- i (4°- (2.6.31) 

=o (2.6.32) 

hsO-e"'")-v!,-vi8 -t- ^=^<1^5 

(2.6.33) 

+ ^ <■" -,?7 -^.] = PR [- '"p“i 

+ =° (2.6.34) 

where the right-hand sides of (2.6.31,33) indicate the 
steady and frequency terms, respectively, of the enclosed 
quantity. The third equation in the set comes from applying 
the boundary condition to at the lower surface of 

the (n+l)st blade 

(2.6.35) 

V>5 (^s = o 



(2.6.36) 



Solving equations (2.6.31) through (2.6.36) for the Sauer 
functions permits determination of the pressure difference 
over the last portion of the blade 
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PL: ^Cp (2.6.39) 



PI; (2.6-.40) 
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2.6.6 Results and Conclusions. Figure 2.6.2 shows that, 
as with the finite-cascade solution, there is considerable 
disagreement among the several existing solutions for small 
interblade phase angles, but even for those angles the 
present solution and Kurosaka's are in agreement. The 
pressure-difference expressions derived above for the infinite 
cascade agree exactly with those of Kurosaka [6] over the 
entire plunging blade and up to ^ 2yj<i vn. for the pitching 
one. Over the last portion of the pitching blade the two 
solutions are slightly different analytically but are in 
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Figure 2.6.2 Comparison of the inf inite -cascade 
■ solution with existing solutions. 



numerical agreement as seen in figures 2.6.3 and 2.6.4; 
the analytical discrepancy could not be resolved. All of 
the solutions are in accord for the larger phase angles. 
Continuing with the conclusions begun in section 2.5.6, it 
can be added that: 

4. For interblade phase angles larger than 15° the 
present method of analysis yields closed-form pressure- 
difference expressions which are in agreement with existing 
solutions for the slowly-oscillating cascade. 

2 . 7 Limiting Cases of the Finite and Infinite-Cascade 
Solutions . 

2.7.1 The Steady Cascade. All terms of the plunge solu- 
tions are frequency terms; therefore if k is set to zero, 
the pressure-difference distribution goes to zero for the 
entire blade, which is, of course, what should happen for a 
stationary flat plate aligned with the flow. 

If k is set to zero in the pitching case, then one 
obtains for the isolated airfoil 




and so both the finite (2.5.27) and infinite (2.6.17) 



solutions becomes 
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Figure 2.6.3 Comparison of the infinite-cascade 
solution with existing solutions. 
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wh 62 TS the phase— angle teriu has been ignored because it is 
meaningless in the steady case. Equation (2.7.2) is the 
classical Ackeret solution [27] for a flat plate at an angle 
of attack equal to one radian, which is the maximum up-pitch 
position. 

2*7.2 The Single-Blade "Cascade" . If either cascade is 
"squeezed" down to a single blade by setting the spacing 
parameters m and d to zero, then the solutions (2.5.26,27) 
and (2.6.16,17) reduce to those of the isolated blade 



where the phase term has again been ignored.- 

2.7.3 The Infinitely-Spaced Cascade. The isolated-blade 
solution is again recovered from equations (2.5.27) and 
(2.6.17) as the cascade spacing is allowed to get large. 

If the spacing parameters m and d approach infinity in 



then equation (2.7.3) again results. The significance of 
these two terms is discussed further in the next section. 

2.8 The Unique Inflow Condition. A consideration of the 
infinitely-spaced cascade in section 2.7.3 leads one^ to 
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and 



remain equal. 



^The author is grateful to Professor M. F. Platzer for 
first asking this question. 
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For which cascades does the relationship 



lic: 

(2.7.4) 

lid? AS seen from equation (2.5.27), for such a cascade 
te pressure-difference over the inflow portion of the blade 
i merely that of the isolated pitching blade times the 

he tor i\- 

Let €>< be the Mach angle defined by 






(2.7.5) 



hen it is not difficult to show that for any cascade having 
he stagger angle 



^ \ (2.7.6) 



jquation (2.7.4) will also hold. Said another way, all 

aetDers of an infinite family of pitching cascades with the 

stagger angle given by (2.7.6) - and with no restrictions 

T-i 1 1 hflve exactly the same inflow pressure 
Dn blade spacing — will have exac y 

distribution, i.e. 



PI: 



,C, = -2C?,U- e"'*) 



(2.7.7) 



This is termed, by the author, the unique inflow condition 
for the slowly-oscillating cascade. Figure 2.8.1 depicts 
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Figure 2.8.1 The unique inflow condition. 



three different cascades, each with the same stagger angle, 
satisfying (2.7.6). The pressure distribution for each 
was obtained via Platzer's method— of —characteristics pro- 
gram. The unique inflow condition is evidenced by the 
exact agreement of the pressures for all three cascades 
over the first portion of the blade. 

2.9 Some Comments on Resonance . For a given reference 
blade in an infinite cascade, resonance occurs when aero- 
dynamic disturbance waves generated by upstream blades are 
in phase with the motion of the reference blade [12] . 
Samoylovich [28] shows that at the resonance condition the 
linear cascade experiences zero aerodynamic damping. The 
development that follows is based on Verdon's [12] excellent 
explanation of the resonance phenomenon and provides a means 
of delineating what he calls the subresonant and superresonant 
regions of the yu-te plane. 

Referring to figure 2.9.1, also taken from [12], it is 
seen that the leading edge of each blade creates disturbances 
which propagate in all directions at the freestream speed 
of sound a^Q . These disturbances are swept downstream by 
the freestream flow U , and there are exactly two resultant 
waves which propagate along the cascade leading edge in 
the direction opposite to the direction of rotation and at 
velocities and V 2 . From the law of cosines one has 







(2.9.1) 




Figure 2.9.1 The backward-traveling waves 

causing resonance. 
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Figure 2.9.2 The subresonant and superresonant regions. 
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The times required for a wave from an upstream 

blade (n — -1, -2, ...) to reach any given reference blade 
(n = 0) are 






= [^yn {vn-v2y34lij 



(2.9.2) 



The reference blade lags an upstream blade by 




n 






(2.9.3) 



and resonance occurs when the upstream waves are in phase 
with the reference-blade motion, i.e. when 
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. (2.9:4) 



or when 







|yi(v«-VZyScJ)j 



(2.9.5) 



For a fixed reduced frequency k the superresonant inter- 
blade phase angles lie between and yUj / shown in 

figure 2.9.2, while the subresonant phase angles lie to 
either side. The superresonant region is important because 
Verdon is of the opinion that, in that region, neither a 
finite-cascade solution nor a passage-approach infinite- 
cascade solution will adequately approximate the true 
infinite-cascade solution. In an effort to shed some light 



on the subject, some comparisons are presented in figures 
2.9.3 through 2.9.6. For cascade A, with k = 0.05 , the 
superresonant range of phase angles is approximately -7° 
to -2° , and so -5° was examined. Figures 2.9.3 and 2.9.4 
show that Platzer's finite solution and the author's passage- 
approach infinite solution are in agreement with Verdon's 
solution. For cascade B, however, Verdon's solution is 
significantly different from the other two over the last 
portion of the blade for the superresonant phase angle -3° . 
Figures 2.9.5 and 2.9.6 show some high-frequency comparisons 
for the phase angle -90° , which is in the superresonant 
regions of both cascades A and B. In these two cases the 
two solutions — one finite and one infinite — are in accord. 
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Figure 2.9.3 Low-frequency superresonant 
pressure distributions. 





Figure 2.9.4 Low-frequency superresonant 
pressure distributions. 
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Figure 2.9.5 High-frequency superresonant 
pressure distributions. 
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Figure 2.9.6 High-frequency superresonant 
pressure distributions. 



3. THE NONLINEAR CASCADE 



3.1 General . In two papers Teipel provides the numerical 
tools for the analysis of low supersonic flow past a single 
oscillating airfoil of finite thickness. In one paper [23] 
he presents an unsteady, nonlinear characteristics method 
to determine the continuous variation of flow variables. 

In the other [24] he treats the oscillating shock waves 
present in the nonlinear problem to permit determination 
of the discontinuous jumps in the flow variables across the 
shocks emanating from the single blade. Using these tech- 
niques as well as some computational and theoretical insight 
provided by. Chadwick [5 ] , who very ably translated and applied 
Teipel 's work, the present study is aimed at finding the 
pressures, steady and unsteady, on the second blade of a 
parabolic-arc cascade. The cascade oscillates in pitch 
alone, and the pressures are sought over the entire upper 
blade surface and over the lower surface up to the point of 
impingement of the trailing-edge shock from the preceding 
blade (figure 3.2.2). These surface pressures are computed 
by means of a FORTRAN IV computer program prepared by the 
author for the IBM- 360 computer. The program listing appears 
at the end of this work immediately after Appendix A. What 
follows is an orderly presentation of the theory underlying 
the computer program and the computational techniques used 
in implementing that theory. Wherever possible the text 
notation will coincide with that used in the program. 



3 . 2 Cascade Geometry . 



3«2.1 Blade Shape . The blade shape was chosen to approx- 
imate the blades in the test cascade employed by Fleeter 
and Riff el [25] who attempted to measure steady and unsteady 
pressures in a pitching cascade. The airfoil used in the 
computer model has a flat upper surface and a lower surface 
conforming to the parabolic arc given by 

-X (3.2.1) 

where chordwise distance x and midchord thickness 'C are 
normalized by chordlength. Figure 3.2.1 shows the computer 
model for Tr = 0.03 superimposed on Fleeter 's test blade 
and lists the surface coordinates of the test- blade in per- ' • 
cent chord. The greatest discrepancy between the two shapes 
occurs on the aft portion of the blade in the region behind 
the lower-surface shock reflection; however, the lower- 
surface pressures are not computed in this region. 

3.2.2 Blade Spacing. The blade spacing parameters are 
input variables for the program and are discussed in more 
detail in section 3.5.1. The physical spacing of the blades 
of the computer-model cascade — dubbed Cascade T — is shown 
in figure 3.2.2. 

3 . 3 Teipel's Method of Characteristics. The underlying 
insight of Teipel's treatment of the unsteady problem is 



I 1 1 1— 

y 0.1 0.2 0.3 0.4 

0.00 




Leading edge to midchord. . 



y ♦ 




Midchord to trailing edge. 





UPPER 


SURFACE 








L0\>raR 


SIRFACE 




★ 


* 


ic 






★ 


* 




ir 


X 


y 


X 


y 




X 


y 


X 


y 


0.00 


-0.26 


49.35 


0.00 




0.00 


-0.26 


50.51 


-2.97 


5.21 


-0.14 


55.65 


0.10 




6.29 


-1.08 


56.83 


-2.34 


11.53 


-0.16 


61.95 


0.19 




12.59 


-1.65 


63.15 


-2.60 


17.84 


-0.16 


68.25 


0.29 




18.90 


-2.12 


69.47 


-2.27 


24.14 


-0.15 


74.55 


0.41 




25.21 


-2.49 


75.78 


-1.83 


30.45 


-0.13 


80.88 


0.50 




31.53 


-2.76 


82.06 


-1.37 


36.75 


-0.10 


87.21 


0.46 




37.85 


-2.93 


88.36 


-0.96 


43.05 


-0.04 


93.53 


0.26 




^4.18 


-3.00 


94.66 


-0.62 






99.85 . 


-O.IC 








99.86 


-0.10 




•* 


■ Expressed ir 


1 


percent 


chord 


1 . 




Surface coordinates. 



Figure 3.2.1 Test blade shape and coordinates. 
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i. representation of the velocity potential as the sum 
i a steady component and a small-amplitude -unsteady one 






! depicted in figure 3.3.1. When (3.3.1)' is substituted 
j!:o the transonic unsteady small-disturbance equation 














39C a-t 



4-M 
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(3.3.2) 



id into the unsteady boundary condition 



^For a discussion of the limitations imposed by this 
Juation see Teipel [23] and Chadwick [5]. 
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Figure 3.3.1 The total velocity potential. 




3-t. 



(3.3.3) 
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where y(x,t) describes the oscillating airfoil surface, 
one obtains two separate problems coupled by the steady 
velocity potential 0 ^ . One problem is nonlinear, steady 
and zero-angle-of-attack 



[^V ° (3.3.4) 

(3.3.5) 

9 ^ 



and the other linear and unsteady 






O (3.3.6) 



LV^CX-b) (3.3.7) 

L J 

where 0^^(x,y) is an amplitude function and the factor ©'■ 
is omitted. It is important to note that since the char- 
acteristic directions are determined by the coefficients of 
the highest derivatives, both problems will have the same 
characteristic directions. The state variables^ for the two 



^Tha symbol / a . is used both for the steady state variable 
and the interblade phase angle , and it is normally clear 
from usage which meaning applies . 



P^oblsrns air© ssssntially the steady and unsteady velocity 
components u^, v^, and 

STEADY: yS*- + M* C^-vO Uo (3.3.8) 

Vi (3.3.9) 



UNSTEADY: 



U|, V, 



(3.3.10) 



where 




(3.3.11) 



The characteristic directions for both problems are • 



=± -i- 



(3.3.12) 



where the upper sign refers to the left-running o< charac- 
teristics and the lower to the right-running ^ character 
istics.^ The compatibility relations, however, are not 
the same for the two problems . They are 



STEADY : 






(3.3.13) 



1 2 — 
The symbol yS is used both for the quantity (M -1)^ 

and the right-running characteristic direction. It is 

normally clear from usage which meaning applies. 
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(3.3.14) 
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where 



^o< = 





(3.3.15) 



Equation (3.3.14) requires some additional comment. 

Note first that in addition to the state variables u^^ 
and , 0^ is also present in the unsteady compatibility 

relations. Since 
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■ (3.3.16) 



then 0^ must be found by integration as one moves from 
one gridpoint to another. For the points A, B, and C 
shown in figure 3.3.2, for example 













(3. 


3.17) 


r- 




-V 






(3. 


3.18) 



In performing the actual computations an average of the two 
above expressions is employed. 



The factor in (3.3.14) is straightforward, repre- 

senting the increment in ^ along the yS characteristic. 
At point C in figure 3.3.2 for example one could use the 
approximation 



(3.3.19) 

However is a bit more complex, for while the increment 

dx is along the c< characteristic, the quantity 
represents the change in A with respect to x along the 
crossing yS characteristic at the point of interest. 

Again at point C one could employ 

(Ac-Aa) (3.3’.20) 

3 . 4 Teipel's Treatment of the Os c illating Shock Wave with 
Extension to the Present Cascade . ^ 

3.4.1 The Isentropic Rankine-Hugoniot Pressure Relation . 

Consider a normal shock moving through an undisturbed 

perfect gas at velocity relative to the gas as shown 

in figure 3.4.1, where the hat notation indicates conditions 

behind the shock. The governing continuity, momentum and 

energy equations are 



9 



^Except for u,, v,, 'll, and 'v, , the velocities appearing 
in sections 3.4.1 through 3.4.6 are dimensional, in a brief 
departure from the normalization of all velocities by the 
freestream velocity U . 




Figure 3.4.1 Jump conditions across a moving, 

normal shock . 




( 3 . 4 . 1 ) 



yo C- Vs') = ^ ^ Vy, - Vs - V„') 

P+ = P CVy,-Vs-V„^^ 



(3.4.2) 

(3.4.3) 



Assuming constant specific heats, the speeds of sound are 
introduced, and the enthalpies are eliminated by 



<x = 









K = 



'S-\ 



(3.4.4) 



(3.4.5) 



and, after a good deal of laanipulation, the change in 
the normal velocities across the shock is found in terms 
of the shock velocity and sound speed ahead of the shock 



Vn-V» _ 2 , . Vt 

(X 



a. 



( 3 . 4 . 6 ) 



Equations (3.4.1) and (3.4.6) and the isentropic relation 






P ' <x ^ 



(3.4.7) 



yield the Rankine-Hugoniot pressure relation sought 






(3.4.8) 



Figure 3.4.2 is based on figure 8-3 of [29], and it shows 
that the isentropic approximation introduces no significant 
error for Fleeter 's test cascade or its computer model. 

3.4.2 The Blade-Two Shock Position. The leading-edge 
shock from the pitching blade oscillates out of phase with 
the blade. Let (y) be the shock angle of the steady, 
leading-edge shock as shown in figure 3.4.3. Since the 
shock is curved, varies with distance y from the 

blade, and by convention 

0< "iJo < ^ (3.4.9) 

The shock angle at any point on the oscillating wave can 
be written as a small variation from the steady one^ 

'6 ^ (3.4.10) 

In order to locate the oscillating shock relative to the 
steady one, Teipel introduces the complex displacement 
function 9(y) / also shown in figure 3.4.3, and so the x 

coordinate of a shock point is approximately 

•X- 'Xo COT'tSo + (3.4.11) 



^The symbol gamma used here and in figures 3.4.4 and 
3.4.5 represents the total shock angle. All other uses of 
the unsubscr ipted , unsuperscripted gamma stand for the 
dimensionless ratio of specific heats. 
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Figure 3.4.4 Velocity components and the oscillating 
shock wave. 
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where ere the steady coordinates of the blade's 

leading edge. Taking the total differential of (3.4.11) 
and noting that along the shock 






= — COX ')i 



(3.4.12) 



results in 




(3.4.13) 



thus relating shock angle and shock displacement. 

3.4.3 The Velocity Components Upstream of the Blade-Two 
Shock Wave . Ahead of the shock the velocity components 
have both steady and unsteady parts 



where the barred velocities are dimensional counterparts of 
the dimensionless quantities u^, v^, u^^ and v^ already 
introduced in equations (3.3.8,9,10). However, the flat upper 
surfaces of the blades of cascade T are aligned with the flow 
in the steady position, and so, upstream of the shock 



W.* = U -V VAo u, 



(3.4.14) 



V* ■= Vo X V| 



(3.4.15) 



Uc, = Vo •= o 



(3.4.16) 





I 



Referring to figure 3.4.4 one has for the normal and 
tangential velocity components ahead of the shock 

V„-=TJ SIN ifo + 'U cos Vo v' C u, + Vj cos Vo") ( 3 . 4 . 17 ) 

Vt=XJcos V*-X3 sm Vo v' -► (u, cos V| SiN Vo'j (3.4.18) 

3.4.4 The Blade-Two Shock-Motion-Induced Velocity . 
Relative to the gas ahead of it the shock is moving at 
the velocity 



Vi = -V„ + v' (3.4.19) 

where V is a normal velocity induced by the motion of 
the shock, and is seen from figure 3.4.5 to be 



V'r:^ 

d-t 






(3.4.20) 



or, from (3.4.11) 



V‘ = tk 



S»N Vo ^ 



(3.4.21) 



3.4.5 The Velocity Jump Conditions Across the Blade-Two 
Shock . Using equation (3.4.13) to eliminate V from the 
expressions for the upstream velocity components (3.4.17,18), 
they become 
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Figure 3.4.5 The shock-motion-induced velocity. 



V„=TJ SIN Vo (u -*• V, cos Vo') (3.4.22) 



Vk-UcosVo+C-^siN^Vo~ -t' u, COSVo-V* (3.4.23) 



Note that each of the above velocities has three component 
parts. The first term in each expression arises from the 
steady shock position, the second arises from the oscillatory 
motion of the shock about that steady position, and the 
third and fourth terms arise from the unsteady disturbances 
ahead of the shock which are generated by the first blade's 



motion. 
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The components behind the shock are found by noting 
that the tangential component does not change 

= (3.4.24) 

and by using (3.4.6,19,21,22) to obtain 
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(3.4.25) 



where 



Mn-=M S\N Vo 



(3.4.26) 



It is convenient to once again return to the dimensionless 
x-y velocity components u^ and , obtaining from the 
unsteady terms of (3.4.24,25) 






'U, = -V* V. y>i2. ^ vr\ju, -4- vn.^V| 



(3.4.27) 



V, = -»■ -t Vlj u, -V n.^ V| 



(3.4.28) 
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where the factor e 



has been dropped and 



ikt 
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wij= cos*- SlN*-'8’o 
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h» = — -r^— C COS Z-tfo -»- SIN* 

O'*" • Mn 



Y>2.^ — — tOT^o SlN*-y, 

V+‘ Mn 



!fl3=: Yn^ 



^ U- ^ -^1 



(3.4.29) 

(3.4.30) 

(3.4.31) 

(3.4.32) 

(3.4.33) 

(3.4.34) 

(3.4.35) 

(3.4.36) 



Note that the above coefficients are functions of the 
local shock angle ^ and so are not constant along the 
shock. 

If no unsteady disturbances exist upstream of the shock 
(u^ = = 0) then (3.4.27,28) reduce to Teipel's single- 

blade results. If, in addition, the shock does not oscillate 
(g = dg/dy =0) , then and are identically zero 

indicating that there are no unsteady disturbances behind 
the shock. Equations (3.4.27,28) represent an extension of 
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Teipel ' s work from the single blade to the cascade, and 
they were first obtained by Chadwick for a cascade of 
wedges. It is noted for completeness that the above 
equations would apply equally well to an upper-surface 
shock wave if the signs of m^, m^ , 1 X 2 and n^ were 
reversed. 

3.4.6 The Pressure Jump Conditions Across the Blade-Two 
Shock. An additional extension of Teipel 's work was accom- 
plished by the author in order to develop a more accurate 
method of determining the flow-variable discontinuities 
than that employed by either Chadwick or Teipel. The 
improved method is treated in more detail in section 3.8, 
but it is based on the extension of the Rahkine-Hugoniot 
pressure relation (3.4.8) to the cascade via 'equations 
(3.4.19,20,21), resulting in 




(3.4.37) 



where 




(3.4.38) 




(3.4.39) 




(3.4.40) 
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(3.4.41) 




Mn cot )So 



Equation (3.4.37) relates the pressure jump across the 
shock to the shock displacement and the upstream unsteady 
velocities. In the single-blade case u^ and are 

zero, and Teipel's result is recovered. 

3.4.7 The Boundary Conditions for the Shock Equations . 
In section 3.8, the flow variables behind the shock are 
obtained by the simultaneous solution of the difference 
equations resulting from (3.3.14) and the shock equations 
(3.4.27,28,37). In order to start the solution process, 
then, several quantities must initially be known at the 
point of intersedtion of the shock and blade, i.e. at 
(x^,y^) . From figure 3.4.3 it is clear that 



Also, immediately behind the shock the flow tangency 
condition for the second blade must hold 






(3.4.42) 



- [l -i' A-t (3.4.43) 



With v^ known equation (3.4.28) provides 




(3.4.44) 
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and (3.4.44) used with (3.4.27,37) gives the remaining 
boundary conditions 










<x.-t ^.TCo, Wo) (3.4.45) 




To.* Vl» 



. )v»i-F{nA 
z, 



(3.4.46) 



tt-t, 



Of course for equations (3.4.43) through (3.4.46) the hatted 



Equation (3.4.45) was derived by Chadwick while (3.4.46) 



3.4.8 The Velocity Potential Jump Conditions Across the 
Blade-Two Shock. Landahl [26] shows that the total velocity 
potential (3.3.1) is continuous across the leading-edge 
shock, but this does not imply that either the steady or 
unsteady component is individually continuous. Teipel and 
Chadwick choose to assume that, at the first shock point 



The author, on the other hand, chooses to find the initial 
velocity potential jump at the point ^ 

perfect gas with constant specific heats, it can be shown that 



quantities are just behind the shock at (x^,y^) . 



is the author's result. 




(3.4.47) 
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(3.4.48) 



^ I . 



^ y 



and since Uj^ and 'P/P» are known at tx^,y^) from 
equations (3.4.45,46), then so is . Figure 3.4.6 is a 

plot of computet data for cascade T, and it shows the 
discontinuity in along the shook as a function of 

distance from the blade. Note that the initial discon- 
tinuity is small justifying Teipel’s assumption (3.4.47), 
but equation (3.4.48) eliminates having to make that 

assumption at all. 



3.5 con.trnctina the ct.r.cteristic Grid. Because the 
characteristic directions are the same for both the steady 
and unsteady problem, the characteristic grid can be con 
strutted and employed to find both the steady and unsteady 
velocities. The grids above blade one (zone A) and above 
blade two (zone C) are composed of two families of straight, 
parallel characteristics. The mesh below blade two (zone B) , 
however, is composed of one straight, divergent family (;) 
and one curved family (®0 • 

3 , 5.1 mt.. Straight characteristic Grids Above Blades On e 

and TWO . Since the upper blade surfaces generate no steady 

” It iri -I that, in zones A and C 

. 1 4 - T c; seen frciu ^:nau, 

disturbances, it is seen 






(3.5.1) 
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Figure 3.4.6 The discontinuous jump in the unsteady 
velocity potential across the leading- 
edge shock.' 
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and so the characteristics in those zones are simply Mach 
waves (figure 3.5.1). Each blade has 21 equispaced upper- 
surface gridpoints , and the program requires that the leading 
edge of the second blade be positioned at one of the straight- 
mesh gridpoints. This is accomplished via the input parameters 
K and N which fix the cascade blade spacing according to 

VB = 0.02.5 (3.5.2) 

r 

X& = vs yS -+- O.os- M (3.5.3) 

For cascade T, K and N are 16 and 7 respectively. 

Figure 3.5.1 also shows the characteristic numbering conven- 
tion employed in the program, so that a gridpoint may be 
located by two characteristic coordinates 






(3.5.4) 



Note that the straight mesh is computed as if both blades 
were flat plates. The x-y coordinates of the straight 
mesh are referred to as XI, Yl in the computer program. 

3,5,2 The Curved Characteristic Grid Below Blade Two . 

In zone B (figure 3.5.2) the characteristic grid is non- 
linear and must be constructed iteratively. Figure 3.5.2 
reveals the numbering convention employed for the curved 
mesh. The curved-mesh x-y coordinates are referred to as 
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Figure 3.5.1 The straight characteristic grid. 
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X2,Y2 in the program to distinguish them from those of 
the straight mesh. Before proceeding with the grid con- 
sturction, some useful relations are obtained. 

For an o< characteristic passing through the shock 
from zone A to zone B (figure 3.5.3) 



/ j 



(3.5.5) 



and since 




is zero 






(3.5.6) 



Now consider two points on the same ^ characteristic, 

also illustrated in figure 3.5.3. From the compatibility 

relation 



?v (X,T) = (3.5.7) 

9vCX^K)^^-yuCX,K) 7^^* (3.5.8) 

and from the yB compatibility relation 

^ (X, -f jx (X,<J) = -V* ^ CX^K) (3.5.9) 



and so it must be that 



>CX,T)= ACI^K) 



r 



iT,Z) = yx(.X,K) 



(3.5.10) 
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family is composed 



which proves, in addition, that the 
of straight lines. 

The steady shock angle at any point along the shock is 
found by averaging the angles of the characteristic lines 
before and after the shock, i.e. 



+TANi’'(V''*)] (3.5.11) 

Taking the tangent of (3.5.11) and applying the sum and 
double-angle trigonometric formulae leads to a quadratic 
equation in tan 

TAN =0 (3.5.12) 

The solution of (3.5.12) is rather cumbersome and is 
adquately approximated by the simpler expression 



T^N V.= 



(3.5.13) 



In the computational implementation of (3.5.13), referring 
to figure 3.5.3 again, the characteristic angles at adjacent 
points behind the shock are first averaged together, and 
the result is averaged with the zone-A angle, i.e. 



TAN = 






(3.5.14) 



iHi 



m 




3. 5. 2.1 Computing the Point (1,1). The X2-Y2 coord- 



inates of (1,1) are of course those of the blade-two 
leading edge. yU(l,l) is found from the steady boundary 
condition (3.3.5) and equation (3.3.9) while 7v(l,l) follows 
from (3.5.6). 

3 . 5 . 2 . 2 Computing the Points (1,2) and (2,2) . The 
points (1,2) and (2,2) are established iteratively as 
follov;s. Program subroutines employed are indicated in 
parentheses, and figure 3.5.4 is provided for reference. 

1. Establish the reference point XA on the blade, 
close to the leading edge 

XA= XZC\,0 -V- O.OOS 

2. Compute ^ at XA , which is equal to yM(l,2) , 
using the steady boundary condition (3.3.5) and (3.3.9). 

3. Compute X at XA , which is equal to 
using equation (3.5.6). 

4. (Subroutine SOLVE) Locate (1,2) 'at the intersection 
of a line through (XA, Y2(l,l)) with slope 



-I 

and a line through (X2(l,l), Y2(1,D) with slope 



A 
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Computing subsequent shock-front points (1,1). 



Figure 3.5.4 Computing the shock-front points. 
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5. (Subroutine ITER) Guess at the position of (2,2) 
on the blade 

JZCZ,!^ ■= XZ(\,0 

6. (Subroutine ITER) Increment that guess 

^( 2 , 2 )= ^C2,Z') 0 . 00 » 

7. (Subroutine ITER) Compute JU{2,2) from (3. 3. 5, 9). 

8. (Subroutine ITER) Compute 7^(2, 2) from (3.5.6). 

9. (Subroutine ITER) Solve for X2(2,2) from the 

slope of the characteristic through (1,2) and (2,2) 

X2(2,2)~X^C\,^) _ Z 

-V- X(.2.,2)*'‘ 

10. (Subroutine ITER) Does X2(2,2) = X2(2,2)? 

YES; GO TO 11; NO: GO TO 6 . 

11. STOP. 

3. 5. 2. 3 Computing Subsequent Shock Front Points . 
Subsequent points (1,1) along the shock are also found 
iteratively. The technique follows, and figure 3.5.4 
applies . 

1. Relocate XA on the blade 



XA.-= XAv + bX/\ 



where DXA is based on the first position of XA 



DXA = 



XZ (Z,Z ) - XA 
15 



2. (Subroutine SIIOKPT) Guess Ad, I) 

= AOjX-0 



3. (Subroutine SHOKPT) Increment that guess. 



AU,X) •= XO^X^ -4* 0.001 



4. (Subroutine SHOKPT) Compute ytX(l,I) using 
A at (2,1) which are the same respectively as 
and A (2, 2) 









(2,z) 



5. (Subroutine SOLVE) Locate (1,1) at the intersection 
of a line through (XA, Y2(l,l)) with slope 



xl 

and a line through (X2 (1,1-1), Y2 (1,1-1)) with slope 

2± 
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from the 



6. (Subroutine SHOKPT) Compute yU(l,I) 
steady boundary condition (3.3.5) at XA and equation 
(3.3.9) . 

7. (Subroutine SHOKPT) Does yU(l,I) =p{l,l)l 

YES: GO TO 8; NO: GO TO 3. ' 

8. STOP. 

3 . 5 . 2 . 4 Computing the Curved-Mesh Interior Points . A 
point (I/J) interior to the curved mesh is located (figure 
3.5.5) at the intersection of a line passing through 
(I-1,J) with slope 



I 

and a line through (I /I) with slope 



-I 

X(.X,X)'/2 

3 . 5 . 2 . 5 Computing Subsequent Blade-Surface Points . 
Subsequent points of the form (1,1) on the lower surface 
are computed in a manner exactly analagous to the computation 
of the point (2,2) as described in steps 5 through 11 of 
section 3. 5. 2. 2. 

3.6 Computing the Steady Velocities and Pressures. With 
yU and Tv. known throughout the grid, the steady velocities 
follow directly from (3. 3. 8, 9). The steady pressures on 




Figure 3.5.5 Computing the curved -mesh 
interior points. 



the upper surface of blade two are zero and on the lower 
surface are given at the surface points (I /I) by 



^ - a Uo (3.6.1) 

3.7 Computing the Unsteady Flow Variables at the 
Straiaht-Mesh Grid Points. 

3.7.1 The Blade-One Upper-Surface Grid Points, . In the 
straight mesh the unsteady velocities and velocity potential 
are designated by 



u^^ = UllR + i UllI = VllR + i VllI 



^ - FlR + i 



FI I 



(3.7.1) 



The unsteady boundary condition for blade one (3.3.7) 
allows determination of the velocities 

VllR(I,I) : -1 VllI(I,I) = -k[xl(I,I) - b] (3.7.2) 

Using the ^ compatibility relation (3.3.14) in difference 
form results in a system of linear eguations for the 
unknowns U11R(I,I), U11I(I,I) . Figure 3.7.1 A applies, 

and the equations appear in matrix form in section A.l of 
the appendix. The velocity potential is found by integration 
of the velocities, using equation (3.3.18); see section A. 2 
of appendix A. 

3.7.2 The Straight-Mesh Field Points . As seen in figure 
3.7.1 B, grid points not on the blade require’ use of both 
compatibility relations. The equations resulting for the 
unknown u^ and v^ velocities at each field point (I,J) 
are given in A. 3 of the appendix. Equations (3.3.17,18) 
are averaged to find the velocity potential; see section 

A. 4 of the appendix. 

3.7.3 The Grid Points Immediately Upstream of the Shock. 
Figure 3.7.1 C shows that the shock points lie in the straight 
mesh and therefore the flow variables at each shock point — 
upstream of the shock — are found by means of an interpolation 
algorithm. The distances RI in figure 3.7.1 C are computed 
and any desired quantity Q(1,K) is found from the values 

of Q at the four surrounding points 





Ui,Vi unkno\^m 
^ / 


\ ^hown 


/. \/^ 
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Figure 3.7.1 Computing the straight-mesh and 
curved-mesh grid points. 
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( 3 . 7 . 3 ) 



QU K) ■= Q(X-t,7)->-RR^ Q(X,~7WK^4- QCT,:y-0 

RD 



where 



RRTS'^'n'R'X = ^ RR-J (3.7.4) 

3 . 8 Determining the Jumps in the Flow Variables Across the 
Blade-Two Shock. It is in the determination of the jump 
conditions that the author employs a technique — alluded to 
in section 3.4.8 — that is different from the method used 
by Teipel in the single-blade case. The essence of Teipel's 
approach is to; • . 

1. Assume that the velocity potential does not change 
across the shock at the first point (3.4.47). 

2. Use the shock velocity equations (3.4.27,28), with 
boundary conditions (3.4.42,43,45), in conjunction with the 




(behind the shock) and g at each shock point. 

3. Use equation (3.3.18) behind the shock to recover 
0^ at each shock point. 

In contrast, the method employed in the present study is to: 
1. Use (3.4.46,48) to determine behind the shock 

at the first shock point. 
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2 . 



Use the shock velocity equations (3.4.27,28), the 
shock pressure equation (3.4.37) , the yS compatibility 
relation (3.3.14) and the boundary conditions (3.4.42,43, 
45,46) to solve for “u^^, g and P/P at each shock point. 

3. Use equation (3.4.48) to find 0^ from the pressure 
jump at each shock point. 

The latter method results in a larger system of equations 
at each point, but no assumptions concerning the continuity 
of the velocity potential are required. 



3.8.1 Details of the Jump-Condition Computations . 

Step 2 above requires some additional comment. The shock 
equations (3.4.27,28,37) are applied in difference form from 
one shock point to the next. They are the, last six equations 
in the matrix equation given in appendix section A. 5. The 
relation supplies the first two equations in that matrix 
equation. Note from figure 3.7.1 D that using the yS relation 
at shock point (1,1) requires two tacit assumptions. First, 
the flow variables at the point where the yg characteristic 
through (1,1) leaves the blade are essentially those at 
the point (1,1) immediately behind the shock. Second, the 
flow variables on the y9 characteristic through (1,1+1) 
and just behind (1,1) are about the same as those at 
(1,1) behind the shock. 

In the computer program the flow variables at shock 
point (1,1) , ahead of the shock, have the form 



ui = U12(1,I) - i U12I(1,I) 



(3.8.1) 



while, at the same point immediately behind the shock, 
they are designated 



= URJ(I) + i UIJ(I) 



(3.8.2) 



until stored in the form of (3.8.1). 

3 . 9 Computing the Unsteady Flow Variables at the Curved- 
Mesh Grid Points. 

3.9.1 The Blade-Two Lower-Surface Grid Points. The 

unsteady boundary condition for the second blade provides 
the lower-surface velocities (figure 3.9.1 A). The 

equations are given in section A. 6 of the appendix. The 

Uj^ velocities are found via the compatibility relation 

(section A. 7) and the velocity potential follows from 
equation (3.3.17) (appendix section A. 8). 

3.9.2 The Curved-Mesh Field Points. Both compatibility 
relations are used at a field point, as seen from figure 
3.9.1 B. The resultant matrix equation appears in appendix 
section A. 9. The velocity potential is found from equations 
(3.3.17,18) (appendix section A. 10). 

With the curved mesh completely com.puted the unsteady 
pressures at blade— two, lower— surface grid points are given 
by 




(3.9.1) 
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3.10 Results and Conclusions. 

3.10.1 The Isolated-Blade Case. If the unsteady boundary 
condition along the first blade is set to zero, the second 
blade sees freestream conditions and its unsteady pressures 
can be compared with Teipel's isolated-blade results [23]. 
This is. accomplished in figure 3.10.1, where it is noted 
that a small discrepancy between the two solutions exists. 
This discrepancy is attributable in part to the different 
techniques used by Teipel and the author in solving the 
shock equations (see section 3.8). But it is also noted 
that Chadwick [5] too encountered a discrepancy when he 
applied Teipel's method to the single,, biconvex airfoil. 
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Figure 3.10.2 High-k comparison of the linear 

and non-linear solutions. 

of Fleeter [25] and the numerical results of the computer 
model is good at midchord and fair at the leading edge. 

The program results compare favorably with theoretical 
predictions from second-order Prandtl-Busemann theory and 
a third-order theory due to Ferri [30] which takes entropy 
losses into account. This latter comparison again demon- 
strates that the computer solution is not limited by its 
use of the isentropic Rankine-Hugpniot pressure relation 
(3.4.8). 

3.10.3 The Unsteady Pressure Results. Figures 3.10.4 
through 3.10.12 show the computer results of the present 
study, the experimental data from Fleeter 's investigation, 
and the linear results of Platzer's method-of-characteristics 
program for nine interblade phase angles. The pressure 
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Figure 3.10.3 Blade-two lower-surface 

steady pressures. 



data are presented as an amplitude and phase angle relative 
to the motion of the second blade. The phase angle is 
measured from the point where the second blade attains 
maximum up pitch to the point of peak blade pressure. Along 
the upper surface of blade two the linear and nonlinear 
numerical results agree exactly - because the upper surface 



is flat — providing an additional check of the nonlinear 
program. Along the lower surface the nonlinear effects 
show up as differences between the two predicted pressure 
distributions — differences which vary with interblade phase 
angle. The predicted nonlinear pressure amplitude is, in 
all cases examined, at least as large as the linear predic- 
tions, which is as expected. In figure 3.10.5 { -4.4”) 
the linear and nonlinear amplitude predictions are nearly 
coincident, while in figure 3.10.9 130.6®) they differ 

considerably. Figure 3.10.13 is a plot of the difference 
between the nonlinear and linear amplitude predictions versus 
the turning angles experienced by the flow at the leading 
edge of blade two. It shows that the nonlinear computer 
predictions .dre consistent since the variation from the 
linear solution increases with the flow turning angle. 

That turning angle is given by 

T.Av.= 3* (.1 - COSy^) -V- 6.8® (3.10.1) 

when blade two is in the maximum up-pitch position. Note 
that the large turning angles in figure 3.10.13 result from 
the choice of a one radian maximum pitch amplitude for 
nondimensionalization purposes. 

Comparison of the nonlinear computer predictions with 
Fleeter 's experimental data, however, reveals that the 
theory and the experimental results do not consistently agree. 
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Figure 3.10.4 Blade-two unsteady pressures. 
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Figure 3.10.5 Blade-two unsteady pressures. 
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Figure 3.10.6 Blade-two unsteady pressures. 
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Figure 3.10.7 Blade-two unsteady pressures. 
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Figure 3.10.8 Blade-two unsteady pressures. 
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Figure 3.10.9 Blade-two unsteady pressures. 
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Figure 3.10.10 Blade-two unsteady pressures. 
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Figure 3.10.11 Blade-two unsteady pressures. 
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Figure 3.10.12 Blade-two unsteady pressures. 
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Figure 3.10.13 Nonlinear variations from linear theory. 

In figure 3.10.11, for example, there is adequate upper- 
surface agreement with the theory, but the lower-surface 
pressure-amplitude comparison is poor. In fact, in no case 
is there agreement between data and theory for both phase 
and amplitude along the lower surface. The lower-surface 
experimental data appear to be somewhat inconsistent. They 
sometimes vary widely from one point to the next for the 
same interblade phase angle (figure 3.10.6) and they also 
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vary a great deal for rather small changes in interblade 
phase angle (figures 3.10.10,11,12). in addition, in each 
of the figures there is at least one amplitude data point 
for the nonlinear (lower) surface that is below the corres- 
ponding linear— surf ace data point. This phenomenon is 
contrary to expectation because the thickness surface must 
generate a larger magnitude disturbance than the flat surface. 
Certainly, additional experimental data are required to 
adequately evaluate the theoretical methods. 

3 *10 -4 Conclusions . In summary, it is concluded that the 
extension of Teipel's single-blade methods to the nonlinear 
cascade results in: 

1. Steady pressure predictions which are in accord 
with other nonlinear theories and which are supported, to 
some extent, by experimental results. 

2. Unsteady pressure predictions which exhibit expected 
trends relative to linear theory, but which are not consis- 
tently in agreement with experimental data presently available. 
The significant influence of thickness on the unsteady 
pressures is apparent in both the theory and data. 

3. A numerical approach to the unsteady cascade which 
appears to have the potential to treat the effects of blade 
thickness, blade camber, and steady angle-of-attack loading 
to provide a more realistic model of the fan and compressor 
blades in a turbomachine. 



APPENDIX A 



A«1 Equations for the Unknown u i Velocity at Blade-One 
Surface Points . 

1. References: section 3.7.1; figure 3.7.1 A. 

2. Axixiliary constants: 

Al = 21 cm2(0.025)/7v|^ A2 = k2M2(0.025)/X;v, 

A3 = 0.0125 A2 A4 = AS/t^* 

A5 = 1-0.0125 A2 A6 = A5/tC^ 

3. Matrix equation: 




•Al/2 

A5 





U11R(I,I) 






U11I(I,I) 





A3 U11R( 1-1,1) 

A3 U11I(I-1,I) - 



> Al Uiri(I-l,I)/2 - A4 V11R(I,I)+A6 VllR(I-l,I) + 
- Al UllR(I-l,I)/2 - A4 V11I(I,I) A6 Vlll( I-l , D + 



+ A2 F1R( 1-1,1) 
+ A2 FlI(I-l,I) 



Ao2 Equations for at Blade-One Surface Points . 

1. Reference: section 3.7.1; figure 3.7.1 A. 

2. Equations: 

F1R(I,I) = F1R(I-1,I) +{ullR(I,I) + U11R(I-1,I) - 

- [vllR(I,I) 'Ir VllR( 1-1,1 (0.0125) 

F1I(I,I) = F1I(I-1,I) +{ullI(I,I) U11I(I-1,I) - 

- [yiii(i,i) Viii(i-I,i3 
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A. 3 Equations for the Velocities u i and v t at Straight -Mesh 
Grid Points. 



1. References: section 3.7.2; figure 3.7.1 B. 

2. Auxiliary constants! 



Al = 2kM^(0.025)/>, 

A. 

A3 r 1 - 0.0125 A2 
A5 = A3/7vy» 

3. Matrix equation: 



A2 = h^M^(.0.025')/p^ 

A4 = 1 + 0.0125 A2 

A6 = A4/v'/i 
A 



A3 


-Al/2 


-A6 


Al/2 


A3 


0 


A3 


-Al/2 


A6 


Al/2 


A3 


0 



0 




UllRCl, J) 


A6 




UllI(I,J) 


0 




V11R(I,J) 


A6 




VllI(I,J) 



A4 UllR(I,J-l) -Y 
A4 Ulll(I,J-l) - 
A4 U11R(I-1,J) 

A4 UllI(I-l,J) - 



Al UlllCl, J-1) - A5 VllR(I,J-l) A2 FlR(I,J-l) 

2 

- ^ UllRd, J-1) - A5 V11I(I,J-1) -V A2 FlI(I,J-l) 

2 

AL UllI(I-l,J) -V A5 V11R(I-1, J) A2 FlR(I-l, J) 

2 : 

- ^ U11R(I-1,J) A5 VllI(I-l,J) -4- A2 FIICI-I, J) 

2 — 



A. 4 Equations for at Straight -Mesh Grid Points . 

F1R(I,J) = -|[f1R(I, J-1) F1R(I-1,J^ ^ 0.0125 UllR(I,J) -t” 

^ 0.00625 |ullR(I, J-1) U11R(I-1,J^ -V 

0.00625 [yllR( I, J-1) - V11R(I-1, j3/>v^ 

FlI(I,J) = i[FlI(I, J-1) -t- F1I(I-1 ,j 3 + 0.0125 UllI(I,J) 

-1. 0.00625 |ulll( I, J-1) U11I(I-1, J)] + 

0.00625 [vll I (I, J-1) - VllKl-l, j3 
1. References: section 3.7.2; figure 3.7.1 B. 
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A, 5 Equations for the Jmp Conditions across the Blade-Two 
Leading-Edge Shock . 

1. References: section 3.8,1; figure 3.7,1 D, 



2. Auxiliary constants: 

Ai = 2/ [j!^(i,i) + Aa,i-i)] 

A2 = 4kM2[x2(l,I) - X2(1,I-1)]/ [>(1,1) i->(l,I-l3 
A3 = A2 k/2 

DY = [y2(1,D - Y2(1,I-i3/2 

3. Matrix equations: 



1 


-A2/2 


Al 


0 


-A3/2 


0 


0 


0 




URJ(I) 


A2/2 


1 


0 


Al 


0 


-A3/2 


0 


0 




UIJ(I) 


1 . 


0 


0 


0 


0 


0 


-m]^/DY 


m 2 • 




VRJ(I) 


0 


1 


0 


0 


0 


0 


-m 2 


-m^/DY 




VIJ(I) 


0 


0 


1 


0 


0 


0 


-n^/DY 


^2 




FRJ(I) 


0 


0 


0 


1 


0 


0 


-n2 


-n-^/DY 




FIJ(I) 


-•^M^ 


0 


0 


0 


0 


^M^k 


-Pl/DY 


P2 




GR(I) 


0 


-tM2 


0 


0 


-•JM^k 


0 


■^2 


-p^/DY 




GI(I) 



URJ(I-l) + A2 UIJ(I-l)/2 -► Al VRJ(I-l) + A3 FRJ(I-l)/2 ■ 
UIJ(I-l) - A2 URJ(I-l)/2 Al VIJ(I-l) -t A3 FIJ(I-l)/2 
m3[ul2R(l,I) + U12R(1,I-1)] 1114 [vl2R(l , I) + Vl2R(l,I-l3- 

m3[ul2I(l,I) U12I(1,I-1)] m4[vl2I(l,I) Vl2I(l,I-l)] - 

m4|ul2R(l,I) + U12R(1,I-D] U4 [vl2R(l,I) Vl2R(l , I-l)] - 

m4[ul2I(l,I) + U12I(1,I-D] ^n4[vl2I(l,I) ->• Vl2I (1 , 1-l )] - 
(P 3 - m2)[u12R(1,I) 4 U12R(1,I-1)] + P 4 [vl2R(l , I ) 4 
(P 3 - m2)[u12I(1,I) U12I (1,1-1)] 4 P4[vl2I(l,I) 4 
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A. 5 ( Continued ^ ; 

3, Matrix equation (right hand side continued): 



* 



- URJ(I-l) - GR(I-1 )/dY - m 2 GI(I-l) 

- UIJ(I-l) - mj_ GI(I-1 )/dY m 2 GR(I-l) 

- VRJ(I-l) - nj_ GR(I-1 )/dY - n2 GI(I-l) 

- VIJ(I-l) - nj^ GKI-D/dY + n2 GR(I-l) 

+ V12RCI-1)] 1c>{m2[f 2I(I) + F2I(I-1)^ M^URJCl-l) - 

^ V12I(I-1^ - k^M^[F2R(I) + F2R(I-l)]-i.'3f M^UIJ(I-l) -*• 



★ 

•A- 



-'SM^k FIJ(I-l) - Pi GR(I-1)/dY - P 2 GI(I-l) 
FRJ(I-L) - GI(I-1)/DY -»- p^ GR(I-l) 



A. 6 Boundary-Condition Equations for the v -|^ Velocity at the 
Blade-Two Surface Points. 

1. References: section 3,9.1; figure 3.9,1 A. 

2. Aiixiliary constants: X = X2(I,I) - X2(l,l) 

3. Equations: 

V12R(I,I) = - COSCyu) - k(X-b )SIN(yu) 

V12I(I,I) = - SIN(yu) - k(X-b)COS(y^) 
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A»7 Equations for the Velocity at Blade-Two Surface 
Points. 



1. References: section 3.9.1; figure 3,9.1 A. 

2. Auxiliary constants (also used in A. 9 equations): 

DXA = X2(I,J) - X2(I-1,J) XLA = + ?V(I-1,J^ 

DXB - X2(I,J) - X2(I,J-1) XLB = i[X(I,J) + A(I,J-1)] 

VAR = i[vl2R(I,J) ^ V12R(I-1,J^ DLB = Xl,J) - X(I,J-1) 
VAI = i[vl2I(I,J) V12I(I,J-1)] A1 = 2 k m2 

DVAR = V12R(I,J) - V12R(I-1,J) A2 = k Al/2 

DVAI t V12I(I,J) - V12I(I-1,J) 

SLA r i ?\(I-1,J)'3 SLB = X(I,J-l)*ll 

DLXB = EX(I,J+-1) - X(I, J)]/[X2(I,J-Vl) - X2(I,J)] 

3. Matrix equation: 

1 ^ DLXB DXA _ A2 DXA^ _ Al DXA 
4 XLA 2 XLA 2 XLA 

Al DXA ^ DLXB DXA A2 DXA^ 

2 XLA 4 XLA " 2 XLA 



U12R(I,J) 

U12I(I,J) 



U12R(I-1,J)[1 
U12KI-1, J) fl 



DLXB DXA . A2 DXA^ l 
4 XLA 2 XLA J 



U12KI-1, J) 



Al DXA . 
2 XLA 



_ DLXB DXA ^ A2 DXA^I 
4 XLA 2 XLA J 



U12R(I-1, J) 



Al DX4 . 
2 XLA 



DVAR 




A2 VAR DXa2 




A2 DXA 


SLA 




XLA SLA 


F2R(I-1, J) 


XLA 


DVAI 




A2 VAI DXA^ 


^ F2I(I-1,J) 


A2 DXA 


SLA 




XLA SLA 


XLA 



4. Note: at surface points I = J. 
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A. 8 Equations for at Blade~Two Surface Points. 

L. References: section 3.9,1; figure 3.9,1 A. 

2. AtJCJciliary constants: as defined in A. 7. 

3. Equations: 

F2R(I,J) z F2R(I-1,J) 4* |•i•[ul2R(I,J) 4* U12R(I-1,j3 

^ var/sla^ dxa 

F2I(I,J) = F2I(I-1,J) 4- {itul2I(I,J) + U12I(I-1 ,j3 

4- VAI/SLa} dxa 

4. Note: at surface points I = J. 



A. 9 Equations for the Velocities u i and at Curved -Mesh 
Grid Points . 

1, References: section 3.9.2; figure 3.9,1 B, 



2. Axaxiliary constants: as defined in A.7; 



3, Matrix equation: 



1 ^ DLXB DXA _ A2 DXA^ _ Al DXA 

4 XLA 2 XLA 2 XLA 



Al DXA -f . DLXB DXA _ A2 DXA^ 
2 XLA 4 :CLA 2 XLA 



I .■ DLB _ A2 DXB^ ^ Al DXB 1 

4 XLB 2 XLB 2 XLB SLB 



Al DXB T . DLB __ A2 DXB^ 
2 XLB 4 XLB 2 XLB 



. 1 _ A2 DXA^ 
SLA 2 XLA SLA 



1 ■_ A2 DXA^ 

SLA 2 XLA SLA 



^ A2 DXB^ 

2 XLB SLB 



1 , A2 DXB^ 

SLB ^ 2 :CLB SLB 



(Note: continued on next page.) 
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A. 9 (Continued): 



3, Matrix equation (unknown vector and right -hand- side 
matrix continued) : 



U12R(I,J) 




U12R(I-1,J) [ 


U12KI, J) 




Ul2I(I-l,J) [ 


V12R(I, J) 




U12R(I,J-1) [ 


V12I(I,J) 




U12I(I,J-1) 1 



DLXB DXA . A2 DXa2"1 


* 


4 XLA 2 XLA J 


1 ^ 


DLXB DXA . A2 DXA^I 


1 « 


4 ^ 2 XLA J 


1 


DLB . A2 DXB^I 


1 + 


4 XLB ^ 2 XLB J 




DLB ^ A2 DXb2] 


1 - 


4 XLB 2 >CLB j 


1 



U12KI-1, •+ Vl2R(I-l,J) 

2 XLA 

- U12R(I-1,J) A?- .■^-^'V12I(I-1,J) 

-mi2I(I,J-l) -V V12R(I, J-l) 

2 XIJ3 

— U12R(I, V12I(I, J-1) 

^ ’ 2 XLB 



f- ^ 


. A2 DXA2 1 . 


L SLA 


2 XLA SLaJ 




J. A2 DXa 2 n ^ 


L SLA 


2 ]a.A SLAj 


r+— 


_ A2 DXb 2 1 a. 


L SLB 


2 XLB SLbJ 


fi-JL 


_ A2 DXb2 1 ^ 


L SLB 


2 XLB SLbJ 



F2R(I-1,J) 
^ F2I(I-1,J) 
-h F2R(I,J-1) 
4- F2I(I, J-1) 



A2 DXA 
XLA 

A2 DXA 
XLA 

A2 DXB 
XLB 

A2 DXB 
XLB 
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A«10 Equations for at Curved-Mesh Grid Points. 

1. References: section 3,9.2; figure 3,9.1 B. 

2. Auxiliary constants: as defined in A. 7. 

3. Equations: 

F2R(I,J) = i[F2R(I-l,J) + F2R(I,J- 

4- ^^U12R(I,J) U12R(I-1,J) J}}dXA 

+ ^[U12R(I,J) 4- U12R(I, J-1) - V12R( I , ^ 

F2I(I,J) = i|P2I(I-l,J) 4- F2I(I, J-1)] 

4 ^[ul2I(I,J) -4 U12I(I-1,J) 4> V12I(I,J)^-^V12 I(I-1,J) ^DXA 

4 i(ul2I(I,J) 4- U12I(I,J-1) - V12KI, J) V12KI, J-l)|nXB 
t SliB 
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